Numerical modelling of blast-induced crack propagation in pre-fractured rock

This paper presents a series of numerical simulations to investigate blast-induced crack propagation in rock specimens embedded with pre-existing fracture. The simulation is conducted using the particle-based discrete element method, which can effectively capture many important geomechanical processes such as solid vibration, growth and fracture sliding. Particles are assembled based on grain size distribution Beishan granite Voronoi tessellation algorithm, parameters particle contacts calibrated against laboratory testing results. represented model smooth-joint contacts. In this study, three sets models established explore effects stress state, length orientation damage evolution specimen subject blast load. first set models, we observe that when horizontal σx not equal vertical σy, maximum slip positively correlated differential stress, i.e. Δσ= − further leads generation different numbers secondary cracks. second shows number fractures grows quadratic power function increased Furthermore, extensive tensile regions emerge around tips fractures. third focuses sensitivity analysis orientation, showing clockwise rotation from direction, decreases significantly but they still concentrated locally at research findings our have implications for assessment excavation damaged zone properties nuclear waste repositories built drilling blasting method.